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matrices of the remaining sensors a so-called bounding sensor is 
calculated, which provides a lower bound to the estimation error 
(see Section [TV}. Both contributions are combined in a branch-and- 
bound algorithm for efficiently traversing the tree for the optimal 
sensor schedule. 
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Abstract — In the considered linear Gaussian sensor scheduling prob- 
lem, only one sensor out of a set of sensors performs a measurement. 
To minimize the estimation error over multiple time steps in a com- 
putationally tractable fashion, the so-called information-based pruning 
algorithm is proposed. It utilizes the information matrices of the sensors 
and the monotonicity of the Riccati equation. This allows ordering sensors 
according to their information contribution and excluding many of them 
from scheduling. Additionally, a tight lower is calculated for branch- 
and-bound search, which further improves the pruning performance. 

Index Terms — Kalman filtering, linear systems, sensor scheduling, 
sensor networks, stochastic optimal control 



I. Introduction 

In sensor systems consisting of multiple sensors, e.g., in sensor 
networks, most information about an observed process is obtained 
if each sensor performs measurements at each time instant. But due 
to constrained communication bandwidth, simultaneously sending all 
acquired data may not be possible. Furthermore, to save energy, 
sensors should be switched off if their measurements are not required. 
To increase the operational lifetime of such sensor systems, the 
measurement rate should be as low as possible, which on the other 
hand increases the estimation error. Sensor scheduling constitutes a 
promising solution to this trade-off. Here, only one sensor per time 
step is activated in order to minimize a function of the estimation 
error over a finite time horizon. 

In this paper, the sensor scheduling problem for discrete-time 
linear Gaussian systems is studied. Determining the optimal sensor- 
schedule, i.e., the sequence of active sensors over the considered 
time horizon that minimizes the estimation error, requires traversing 
a search tree as shown in 1 1 1. The depth and branching factor of the 
tree are equal to the time horizon length and the number of sensors, 
respectively. Each path from the root to a leaf node of the tree corre- 
sponds to one possible sensor schedule. The tree-search corresponds 
to solving a binary integer program and thus, sensor scheduling is 
NP-hard. Especially for long time horizons and/or a large number 
of sensors, an exhaustive tree-search can be avoided by employing 
pruning techniques. Existing pruning techniques can be classified into 
suboptimal and optimal methods. By employing optimal pruning [2|, 
(U, ID, deleting the optimal sensor schedule is impossible, but many 
complete sensor schedules may still be evaluated. Drastic savings in 
computational demand arise from suboptimal pruning [5], |6|, but 
ending up with the optimal schedule is no longer guaranteed. The 
sub-optimal methods in |4|, [7|, |8| utilize convex relaxation and 
optimization instead of traversing the tree. In |9|, |10|, the sensor 
scheduling problem is formulated as partially observable Markov 
decision process and solved approximately by means of discretizing 
the belief space. 

In the following, an optimal pruning technique named information- 
based pruning (IBP) is introduced aiming at pruning complete sub- 
trees as early as possible. To achieve this, two contributions are made: 
First, a partial ordering of sensors based on the so-called sensor 
information matrix (see Section |ill). This allows excluding sensors 
from tree search merely on the basis of the measurement matrix 
and noise covariance matrix. Second, with the sensor information 



II. Problem Formulation 
The dynamics of the observed process are given by the linear model 



k.fc+1 



= A fe • x k + w k 



(1) 



where k — 0,1, . . . ,N — 1 is the discrete time index and N is the 
time horizon lengthrlA finite set S of sensors is considered, where 
measurement zjj. from sensor i S 5 = {1, . . . , S} is related to the 
latent system state x k £ R™" via the linear measurement model 
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Both Afc and H\ are time-variant matrices. The noise terms w k 
and v* h are zero-mean white Gaussian with covariance matrices 
C k and C k ' 1 , respectively, with C k '' being positive definite. A 
measurement value i k of sensor i 6 S is a realization of z_ k . The 
initial system state x ~ Af(x ; x , Cg) at time step k = is known 
and assumed to be Gaussian with mean vector x and covariance ma- 
trix Co . A sensor schedule is indicated by u — (wo, iti, . . . , ujv-i), 
where an element Uk encodes the index of the sensor to be scheduled 
for measurement at time step k, i.e., if sensor i is scheduled at time 
step k then itfe = i. 

For linear Gaussian systems, the Kalman filter is the optimal 
estimator of x k in a minimum mean-square sense. For a given sensor 
schedule u, the state covariance matrix C% evolves according to the 
Riccati equation (partly in information form) 
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where the function n(C%) propagates the state covariance C k to the 
next time step given the measurement of sensor i. Further, MJ, is 
the symmetric and positive semi-definite sensor information matrix. 
It subsumes the information contribution of sensor i at time step k. 
The aim of multi-step sensor scheduling is to determine the sensors 
for the time steps k = 0,1, . . . ,N — 1 such that the uncertainty of 
the state estimate and thus, the estimation error is minimized. For this 
purpose, the optimal sensor schedule u* — (uq, . . . , Ujv-i) results 
from solving 



with 



u* — arg min{ J(u)} 

u 

N 

J(u) = X>(W fc CE(u)W 



(2) 






where W^ is a positive semi-definite weighting matrix. By means 
of Wfc, an adaptation to specific scheduling objectives is possible. 
For instance, with Wfe = I for k — 1,...,N the focus is on 
minimizing the average estimation error, while for W^ = for 
k — 1, . . . , N — 1 and Wjv = I merely the terminal estimation error 
is minimized. The cost function gu can be the trace, determinant, 
or the maximum eigenvalue of the weighted state covariance matrix. 
Thus, gk maps the weighted estimation error or uncertainty subsumed 
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'Random variables are denoted by lowercase bold-face, vectors by under- 
lined, and matrices by uppercase bold-face letters. 
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Fig. 1. Search tree for time horizon length N = 2 and 5 = 2 sensors with 
root node CS. 



in W fc C|(M)Wfc , k = 1, . . . , iV to a scalar value. J(«) is the foto/ 
cosf of the schedule u. 

Solving the optimization problem in |2]( corresponds to finding the 
path with lowest cost from the root to a leaf node in a tree like 
the one depicted in Fig. [T] Obviously, enumerating all possible paths 
is computationally infeasible for many sensors and/or a long time 
horizon. 

III. Partial Ordering of Sensors 

The first step of the proposed IBP is based on the monotonicity 
property of the Riccati equation. 

Theorem I For any i g S, positive semi-definite matrix Wfc, and 
covariance matrices C%, 0% with C% >; C%, i.e., C^ — Of is positive 
semi-definite, it holds: 

(Monotonicity) n{C%) b n{C%) , 

(Cost ordering) S*(C£) > g k {C x k ) , 

(Symmetric weighting) W fc C£W^ X W fc CfcW£ . 

For a proof of the monotonicity see Lemma 2 in II II and for 
the cost ordering property see |12|. The monotonicity implies that 
the positive semi-definite ordering between covariance matrices will 
not change by applying the Riccati equation. This ordering is also 
not affected by multiplying the covariance matrices with the positive 
semi-definite weighting matrices Wfc in J2l. 

Furthermore, the ordering property automatically implies an or- 
dering of the (scalar) cost functions. But pruning nodes merely on 
the basis of the scalar costs gu without considering the ordering of 
the covariance matrices automatically leads to suboptimal pruning. 
This follows from the fact that the cost ordering is only a necessary 
but no sufficient condition for the positive semi-definite ordering 
between C\ and C%. Hence, at each node of the search tree, the 
ordering of covariance matrices needs to be determined explicitly for 
the next time step in order to ensure optimal pruning. On this account, 
the optimal pruning procedure proposed in |2| evaluates the Riccati 
equation for each sensor, which is computationally demanding for 
large state spaces or high dimensional measurements. The proposed 
pruning technique instead allows determining the order of covariance 
matrices without evaluating the Riccati equation multiple times. 

Theorem 2 (Order of Sensor Information Matrices) 

Given the covariance matrix C% and the sensor information matrices 
MJ, and TS/lj, for two sensors i,j £ S such that 



Ml y M{ 



(3) 



PROOF. Adding (C%) x to both sides of J3j and inverting both sides 
yields ((C^ 1 + M l fc ) _1 X ((CJ)" 1 + M£) -1 . Both sides are 
now multiplied from right with Aj and from left with A*,, Adding 
Cf results in n(C%) X rj(C%) . □ 

Thus, for checking the positive semi-definite ordering of covariance 
matrices, merely the ordering {3} of sensor information matrices for 
two sensors i and j needs to be determined. Selecting sensor i will 
then provide a smaller covariance matrix at time step k + 1. Due to 
the monotonicity of the Riccati equation, this even holds for arbitrary 
sensor schedules from time step k + 1 on. 

Corollary 1 (Pruning via Sensor Information Matrices) Suppose 
that the sensor schedule u — (wo,Wi, ...,Uk—i) up to time step 
k — 1 results in the covariance matrix G%(u ) . If for two sensors 
Uu and Uk the order 

Ml" X rvC= 

holds, then Vu — (Wfe.fi, . . . , ujv-i) there exists a schedule u = 
(uk+i, ■ ■ ■ , Ujv-i) such that J (u) > J (u), where u = («', Uk, u') 
and u — (u ,uu,u ) . 

PROOF. Due to Theorem [2] at least the sensor schedule u" = u" 
yields C%(u) >^ C%(u) for all k. The cost ordering according to 
Theorem [I] then leads to J(u) > J(u) . □ 

Without evaluating the Riccati equation at time step k and only 
by comparing the sensor information matrices of sensor ut and 
Uk, it can be decided that only sensor ut needs to be considered 
for determining the optimal sensor schedule, while the covariance 



matrices C 



fe+i 



(«), 



,Cff(u) need not to be determined for any 



thenn{C%) <r 3 (C\ 



sequence u . Hence, the complete sub-tree of sensor Ufe can be 
pruned. 

A. Order of Sensor Information Matrices 

When comparing sensor information matrices (or equivalently 
covariance matrices), the difference is in some cases indefinite, i.e., 
it is not determinable, if one information matrix is "larger" than the 
other. This is due to the fact that the order relation X of positive 
semi-definite matrices leads to partial orders. Thus, in general, not 
all sensors can be pruned at a specific time step. 

Scalar systems, i.e., the dimension of the state x is one, form an 
exception. Here, the partial order becomes a total order. Per time 
step, it is now possible to prune all sensors except of one, which is 
equivalent to selecting the sensor that minimizes the state variance 
at each time step. This greedy strategy automatically leads to the 
optimal sensor schedule. Similar results can be found in 1 13 1 for the 
case of scheduling multiple scalar systems to a single sensor. 

B. Branch-and-Bound 

Even if pruning based on sensor information matrices reduces the 
number of possible sensor schedules, the remaining number of nodes 
in the search tree may still be large due to the partial order. To further 
prune the tree, the proposed technique is combined with branch-and- 
bound search. Branch-and-bound is a common search technique for 
classical decision problems like traveling-salesman or knapsack. The 
basic idea is to assign a lower bound of the achievable total cost value 
to any visited node. By means of these bounds, the tree is further 
expanded (branched), whereas nodes with a smaller lower bound are 
considered more promising to lead to the optimal sensor schedule and 
thus are expanded first. Sub-trees are pruned if their lower bound is 
larger than the total cost value of an already completely evaluated 
sensor schedule. 



Algorithm 1 IBP branch-and-bound algorithm. The algorithm is initialized with J m i n = oo. 



For a given sensor schedule (uo, 
if leaf node, i.e., k = N then 



. ,Uk-i) do: 



Jmin <— J(uq, . . . 

else 

hi ^— child(Ufc„i) 
hi' <— prune(W) 
bound(W) 
U' <- sort(W') 
U" *- prune(W) 
expand(W") 
end if 



,un-i) 



II children of sensor u^-i 

II use ordering of sensor information matrices 

// calculate lower bounds for all sensors in hi' 

II sort in ascending order based on lower bounds 

// use lower bounds and J m i n 

// call Algorithm 1 for remaining sensors Ut, £ hi" 



For a particular node that was reached during the search by 
employing the sensor schedule (uo, . . . ,Wfc-i), the total costs can 
be written as 



J(u) — J(uo,- ■ ■ ,Uk-l)+J(Uh, ■ ■ ■ ,Un-i) 



(4) 



where only the value of the first summand is evaluated and thus 
known. While the value of the second summand is not calculated 
yet, a lower bound can be easily assigned by exploiting the sen- 
sor information matrices as shown in the next section. Based on 
this bound, the sub-tree corresponding to the remaining schedules 
(v,k, ■ ■ ■ ,um-i) can be pruned, if the lower bound is larger than 
a global bound J m m, which is the total cost of the currently best 
completely evaluated sensor schedule. J m j n forms an upper bound 
of the optimal cost. Obviously, the closer the lower bound to the 
true value of J(uk, . . . , wjv-i), the earlier complete sub-trees can 
be pruned and thus, the better the pruning performance. 

IV. Bounding Sensor 

Due to the cumulative structure of the total costs J(u) with non- 
negative summands, a simple lower bound can be obtained by setting 
the second summand in (|4j equal to zero. In the following, this 
simple bound is referred to as zero bound (ZB). This bound can be 
interpreted as expecting a complete reduction of uncertainty for each 
time step k, ... ,7V — 1. Such a reduction can merely be achieved 
by a sensor information matrix with infinite trace or determinant. 
Obviously, such a sensor information matrix provides the coarsest 
lower bound possible. 

A significantly tighter lower bound can be derived by means of 
a so-called bounding sensor with sensor information matrix M^ for 
which the covariance matrices evolve according to 



C n +i — A, 



((ciy 1 



+ M T 



A^ + CI 



(5) 



for n = k, . . . , N - 

step n, this special information matrix fulfills 



For each time 



(6) 



for all sensors u n G 5, where M n is as close to MJ^" as possible. 
According to Theorem [2] for covariance matrices obtained from |5} 
holds Cfj +1 ;< C^ +1 (un) for all time steps n = k, . . . ,N — 1 and 
sensor schedules (ttfe, . . . , mjv-i). 

To prune the tree at an arbitrary time step (or level) k, merely the 
cost for the bounding sensor needs to be evaluated for the time steps 
k, . . . , N — 1 in order to obtain the lower bound. This calculation 
can be done in linear time. In comparison, calculating the costs of 
all possible sensor schedules (v,k, ■ ■ ■ ,Un-i) requires exponential 
computation time. 



The complete IBP search algorithm is summarized in Algorithm [T] 
where as basic structure a depth-first search is applied. Breadth-first 
search can be utilized as well. The determination of the bounding 
sensor with information matrix M&, which is required for calculating 
the lower bounds in line 7, is described in the following. 

A. Bounding Sensor for Two Sensors 

Since sensor information matrices are symmetric and positive semi- 
definite, they can be graphically interpreted as ellipsoids similar to 
the covariance ellipsoids that correspond to covariance matrices. In 
the following, the ellipsoid corresponding to a sensor information 
matrix is referred to as sensor information ellipsoid. Unlike covari- 
ance ellipsoids, sensor information ellipsoids have no distinguished 
position. Hence, it can be assumed that all sensor information 
ellipsoids are centered around the origin. 

Based on this interpretation, determining the bounding sensor 
information matrix M^ corresponds to the determination of the 
covering ellipsoid that contains the ellipsoids of all sensor information 
matrices. This problem is similar to the so-called Lowner ellipsoid 
problem (see e.g. |14|). Thanks to the fact that all sensor information 
ellipsoids have the same center, the Lowner ellipsoid problem can be 
significantly simplified. 

At first, the two sensors case is considered. Let vr- 1 ' , V, 2 ' 6 S be 



(i) 



T (2) 



the two sensors with information matrices M. k , TS/l k . For this case, 
a bounding sensor information matrix with minimum determinant, 
i.e., a minimum volume covering ellipsoid, results from solving 
a generalized eigenvalue problem as summarized in the following 
theorem and as depicted in Fig. [2] 

Theorem 3 (Min. Bounding Sensor Information Matrix) The 

bounding sensor information matrix Mfe with minimum determinant 
that fulfills ^ for two sensor information matrices M. k and M. k 
is given by 

-l 



[Pi! 



M fc = (V 1 J ■ max^M*: 1 ' V, V T M^ ; Vj ■ V - " , (7) 
where max( ■ ) is the element-wise maximum of matrices and V is 



the matrix of generalized eigenvectors of MJ. and M 



(i) 



r(2) 



rd) 



(V 



PROOF. At first, the sensor information matrices M. k and M fc 
have to be diagonalized simultaneously, which requires solving the 
generalized eigenvalue problem 



rd) 



M[.' - AM.^ I = 



(2) I 



(8) 



The solution of {5} is given by the eigenvalues Xi,i £ {1,2, ... ,n x \ 
and the matrix of eigenvectors V. This allows diagonalizing Mi. 

(2) 

and MJ. according to 

V — UK!'; ! I, \ ] /\,, I . 

(9) 



v t m;/ ; v 



r(2) 



diag([Ai, 
V = I . 




(b) " 
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Fig. 2. Illustration of the determination of the minimum bounding sensor 
information matrix for the two sensors' case, (a) Sensor information ellipsoids 



of two sensor information matrices Mi, . Mi, . (b) Result of the simultane 



(i) 



r (2) 



r (l) 



(2) 



ous diagonalization of Mj_ and M^ . The covering ellipsoid (red dashed 
ellipsoid) corresponding to the bounding sensor information matrix results 
from taking the maximum eigenvalue (or longest principal axis, indicated 
by the arrows) for each dimension, (c) Transforming back yields the desired 
minimum covering ellipsoid corresponding to M/,. 
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Fig. 3. (a) Average number of expanded nodes, (b) Average deviation of 
greedy scheduling and convex scheduling from optimal scheduling in terms 
of total costs J. 



A sensor network observes the target. It consists of eight sensors 
with measurement matrices 



Determining the minimum volume covering ellipsoid for diagonal 
matrices is straightforward by taking the maximum of each diagonal 
element of the matrices in {9l according to 



maxfV'M 



(i) 



v,v t m1 2) v) 



(10) 



This corresponds to taking the longest principal axis for each di- 
mension or equivalently to taking the maximum of (1, A») for each 
i e {1,2, ...,714 ( see Fig. El (b)). Multiplying (V T ) _1 from left 
and V -1 from right to |l0j reverses the diagonalization and leads to 
the desired minimum sensor information matrix M^ (see Fig. [21(c)). 

□ 



B. Recursive Calculation for an Arbitrary Number of Sensors 

Determining the bounding sensor with minimum sensor informa- 
tion matrix for an arbitrary number of sensors is computationally ex- 
pensive in general, as numerical optimization is required. However, a 
very tight but not necessarily minimum bounding sensor information 
matrix can be efficiently calculated by employing the solution of the 
two sensors' case recursively in a pairwise fashion. In doing so, the 
complexity for determining the bounding sensor information matrix is 
linear with the number of sensors. 

The recursion commences from the initial solution M). according 



to (7]) for the first two sensors 



,(1) „(2) 



G S. For the remaining 



sensors u^' £ S, with i £ {3, 4, . . . , |<S|}, the recursion 



M 



W 



■ njax- (V T M^ 1) V,V T M< ; i) V) -V -1 (11) 



applies, where V is the matrix of eigenvectors of M^ and M 



r<<> 



The final solution M), for i = \S\ is the desired bounding sensor 
information matrix Mfe that fulfills (6|. 

V. Simulation Example 

The proposed IBP is evaluated by means of a simplified target 
tracking example. The state x k = [xk, Xk, y^, Vk] T °f me observed 
target comprises the two-dimensional position [xk,y k ] T and the 
velocities [xk,y k ] T in x and y direction. The system matrix and 
noise covariance matrix of w k of the dynamics model {T} are 



A = I 



and C w = q ■ 1 1 



L> 



2 A 



(12) 



respectively, where <g> is the Kronecker matrix product. In dT2b, T = 
1 s is the sampling interval and q = 0.02 is the scalar diffusion 
strength. Mean vector and covariance matrix of the initial state x 
are x = [0, 1,0, 1] T and Cg = I, respectively. 



H 1 = [1 0] , 

H 3 = [0 1 0] , 

H 5 =[(H^) T ,(H 3 ) T ] T 

H 7 =[(H^) T ,(H 2 ) T ] T 



H 4 



[0 10 0] 
[0 1] 



H 6 =[(H 2 ) T ,(H 4 ) T ] T 



H 



[(H 3 fc ) T ,(H 4 ) 



4nT]T 



The noise covariance matrices C k are diagonal, where the diagonal 
elements are determined randomly from a uniform distribution on 
the interval [0, 1]. For each k the weighting matrix W^ = I, i.e., 
the objective is to minimize the average estimation error. The cost 
functions gk(-) are set to the trace for each k. Six different pruning 
methods are compared for time horizon lengths A^ = 1, . . . , 8: 

GRE Greedy scheduling, i.e., one-step lookahead. 
CVX Suboptimal scheduling based on convex 

optimization |8|. 

Branch-and-bound with zero bound, see Section 

Pruning method proposed in |2l . 



ZB 

COV 

SIM 

IBP 



IV 



ZB plus the order of sensor information matrices. 
Proposed information-based pruning. 
For each time horizon length, M = 50 Monte Carlo runs are 
performed. 

In Fig. [3] (a) the average number of expanded nodes in the search 
tree is plotted. Compared to the total number of nodes, which is 
E!Lol 5 r G C(\S\ N ) with O(-) being the big-0 in Landau 
notation, optimal pruning significantly reduces the search space. The 
much better performance of SIM compared to ZB shows the benefit 
from utilizing the order of sensor information matrices. Additionally 
employing this order can be considered a pre-selection of candidate 
sensor schedules, while branch-and-bound thins out this candidate 
set. IBP provides a much tighter lower bound than ZB and SIM, 
respectively, which further improves pruning. 

The average deviation of the suboptimal methods GRE and 
CVX from optimal scheduling is depicted in Fig. [3] (b). The av- 
erage deviation for a given time horizon length N is defined as 
it? X^i-i i^N,' ~ Jn) w i m Jn being the optimal total cost of 
simulation run i calculated by means of IBP and ,P N , being the 
total cost of pruning method • G {GRE, CVX}. More significant 
deviations from the optimal scheduling especially in case of GRE are 
expected for example in scenarios where additional sensor constraints 
like energy or communication costs (see Section |VI-B| > have to be 
considered, where sensors are temporarily unavailable, or where the 
dynamics and sensor models are nonlinear (consider for example the 
results in |15|, |16|). Regarding the runtime, IBP is outperformed 
only by GRE, while all optimal scheduling algorithms and CVX are 
up to three orders of magnitude slower than IBP. 



VI. Potential Extensions 

The proposed information-based pruning algorithm can be ex- 
tended in many ways in order to apply it to more general sensor 
scheduling tasks and to further improve the pruning performance. 
Three possible extensions are discussed. 

A. Multiple Sensors per Time Step 

So far, merely one sensor per time step is selected for performing a 
measurement. The extension to multiple sensors requires to schedule 
a subset of sensors S' C 5 with \S'\ = S' per time step. For this 
purpose, each subset out of („,) possible subsets has to be evaluated 
at each time step if no pruning is performed. This additional subset 
selection has to be included into the tree search in order to employ 
IBP for reducing the search effort. 

A straightforward realization is to perform nested branches per 
time step instead of a single branch as in the single sensor case. 
Each nested branch can be considered a tree of depth 5" but with 
a branching factor that decreases with each level. In the first level, 
the branching factor is S, in the second level it is S — 1 since the 
sensor selected in the previous branch must not be considered in the 
current branch. Finally, in level S' the branching factor is reduced 
to S — S' + 1. Each path from the root to a leaf node in this subset 
selection tree corresponds to a sensor subset that can be scheduled 
for a time step. At each leaf node, the subset selection tree for the 
next time step begins. Overall the whole search tree has a depth of 
N ■ S' with varying branching factor. For comparison, the search tree 
for the single sensor case has a depth of N with a constant branching 
factor of S. IBP can now be directly employed as in the single sensor 
case. 

B. Sensor Constraints 

In many sensor scheduling problems, there is an additional budget 
constraint of the form 

B-u<b, (13) 



where u — [Uq\ 



£ {0,1}^ is a binary sensor 



scheduling vector with elements u k £ {0,1} for each time step 
k — 0, . . . , N — 1, Such a budget can for example be a maximum 
number of sensors K 6 IN allowed to perform measurements. Here, 
B = [1, 1, . . . , 1] T is a vector of ones and b = K. A further example 
is an energy budget K < £ IN that is assigned to each sensor i. The 
energy budget Ki is reduced by one after each measurement of sensor 
i. In this case, b = [Ki, Ki, ■ ■ ■ , Ks] T and B is a binary matrix, 
where only the elements i, i + S, . . . , i + (N — 1) • S of row i of B 
are equal to one. 

If such sensor constraints are present, IBP cannot be applied 
directly. This is basically due to pruning based on the partial ordering 
of sensors as this ordering does not take the sensor constraints into 
account. Additionally to considering the sensor information matrices, 
it is also necessary to check whether the budget constraint can still be 
fulfilled if a node of the search tree is pruned. Consider the example 
mentioned above, where a maximum number of sensors K € IN is 
allowed to perform measurements. To map IBP on this constrained 
scheduling problem, one possibility is to include a virtual sensor with 
information matrix M = 0. Selecting this virtual sensor corresponds 
to performing no measurement. No budget costs are assigned to this 
virtual sensor, i.e., the elements in B corresponding to the virtual 
sensor are zero. Of course, the information matrix of the virtual 
sensor is always dominated by all other sensors, but selecting the 
virtual sensor will never violate the constraint K. Thus, by ensuring 
that this sensor is never pruned, IBP can be applied directly. 



In contrast to the partial ordering, calculating the lower bound 
based on the sensor information matrices is always directly possible. 
Certainly, the bound will not be as tight as in the unconstrained case. 
Calculating the bound at a particular node in the search tree requires 
to assume the best case where each sensor is available at all times. 
Thus, all sensor information matrices are incorporated into calculating 
the bounding sensor even if some sensors might not be available in 
future. 

It is worth mentioning that the constraint |T3j itself can be 
considered as some kind of pruning. Avoiding the selection of a 
sensor that will violate the constraint corresponds to pruning the 
whole sub-tree of this particular sensor. 

C. Improving the Bounds 

Scheduling algorithms based on convex optimization (see for 
example |4|, |7|, |8|) are sub-optimal, but they can be employed in an 
optimal branch-and-bound algorithm for calculating lower bounds as 
shown in |4|. In many cases, these bounds are tighter than the bound 
proposed in this paper. However, calculating lower bounds based on 
convex optimization is computationally significantly more complex. 
Thus, both bounds could be combined, where the proposed bound 
is used firstly for pruning in order to reduce the need of calculating 
bounds via convex optimization. 

For IBP merely the total cost of the currently best completely 
evaluated sensor schedule is used for setting the global bound value 
Jmin- By means of tight upper bounds, J m i n can be decreased more 
rapidly, which will significantly improve the pruning performance. 
One way to calculate upper bounds is to use Kalman predictions 
without measurement updates. Better upper bounds can be provided, 
if the sensor information matrices are used for calculating a further 
bounding sensor with information matrix M' fc X M^* similar to |6l. 
Therefore, the recursive algorithm proposed in Section |IV| can be 
employed, where max in <JT0J> and \\\\ is replaced by min . 

VII. Conclusion 

The proposed information-based pruning excludes complete sub- 
trees very early, while preserving the optimal sensor schedule is 
guaranteed. Compared to existing optimal pruning methods, the ef- 
fectiveness of the proposed method relies on exploiting the properties 
of the sensor information matrices. This allows pruning without 
explicitly evaluating the Riccati equation and by calculating a tight 
lower bound to the total costs. IBP can be extended for example to 
schedule multiple sensors per time step or to consider constraints like 
energy budgets. 
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